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ABSTRACT 

Stellar rotation periods measured from single-age populations are critical for investigating how 
stellar angular momentum content evolves over time, how that evolution depends on mass, and how 
rotation influences the stellar dynamo and the magnetically heated chromosphere and corona. We 
report rotation periods for 40 late-K to mid-M stars members of the nearby, rich, intermediate-age 
(~600 Myr) open cluster Praesepe. These rotation periods were derived from ^200 observations 
taken by the Palomar Transient Factory of four cluster fields from 2010 February to May. Our 
measurements indicate that Praesepe's mass-period relation transitions from a well-defined singular 
relation to a more scattered distribution of both fast and slow rotators at ~0.6 Mq. The location 
of this transition is broadly consistent with expectations based on observations of younger clusters 
and the assumption that stellar-spin down is the dominant mechanism influencing angular momentum 
evolution at 600 Myr. However, a comparison to data recently published for the Hyades, assumed 
to be coeval to Praesepe, indicates that the divergence from a singular mass-period relation occurs 
at different characteristic masses, strengthening the finding that Praesepe is the younger of the two 
clusters. We also use previously published relations describing the evolution of rotation periods as a 
function of color and mass to evolve the sample of Praesepe periods in time. Comparing the resulting 
predictions to periods measured in M35 and NGC 2516 (^--^150 Myr) and for kinematically selected 
young and old field star populations suggests that stellar spin-down may progress more slowly than 
described by these relations. 
Subject headings: stars: rotation 
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1. INTRODUCTION 

In a seminal paper, Andrew Skumanich (1972) showed 
that stellar rotation decreases over time such that Vrot oc 
i""'^, as does chroniospheric activity, a proxy for mag- 
netic field strength. This relationship between age, ro- 
tation, and activity has been a cornerstone of stellar 
evolution work over the past 40 years, and has gen- 
erated almost as many questions as applications. An- 
gular momentum loss due to stellar winds is gener- 
ally thought to be responsible for the Skumanich law, 
but the exact dependence of Vrot on age is not en- 
tirely understood, and relies on the assumed stellar mag- 
netic field geometry and degree of core-envelope coupling 
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(jKawaled llQSa iKrishnamurthi et all Il997f ). Further- 
more, later-type, fully convective stars appear to have 
longe r active lifetimes than their early-type brethren 
(e.g., iWest et all I2008D . indicating that they are ca- 
pable of generating significant magnetic fiel ds even in 
the a bsence of a standard solar-type dynamo [jBrownind 
l2008f l. The lack of a comprehensive theoretical under- 
standing of the age-rotation-activity relation has not 
prevented the development and use of gyrochronol- 
ogy, however, which is used to determine the ages 
of field st ars b ased on a presumed ag e-rotation rela- 
tion (e.g ■ iBarn cs 2007: Mamaick & Hill enbrand! 120081 : 
iCollier Cameron et al.ll2009l : iBarnesil201Ql , nor of empir- 
ical age-activity relations, which do not always find ac- 
tivity decaying q uite as simply as predicted by the Sku- 
manich law (e.g..lFeigelson et al.ll2004l : iPace fc Pasauinil 
[200l iGiampapa et al.ll200ir 

Mapping out the dependence of stellar rotation and 
activity on age requires the study of stars ranging in 
both mass and age. Statistical constraints on the age- 
rotation-activity r^lation_can_be_d£rivBd_^tm 
field stars fe.g.. iFeigelson et al.ll200l iCovev et~alll2008l : 
llrwin et al.ll2011l ). but the homogeneous, coeval popula- 
tions in open clusters provide an ideal environment for 
studying time-dependent stellar properties. Ideally, rota- 
tion periods for large numbers of cluster members could 
be measured directly from modulations in these stars' 
light-curves due to the presence of star spots. There are 
relatively few nearby open clusters, however, and fewer 
still have had the high quality photometric data needed 
to characterize their members' rotation in this manner — 
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in part because of the sheer difhculty involved in system- 
atically monitoring a large number o f stars over s evera l 
months or more. Studies like that of lSkuman ich' (1972) 
relied instead on measurements of the rotational Dopplcr 
broadening of spectral lines, a technique that has the 
advantage of needing only one observation. Translating 
the resulting Vrot sin i measurements into Prot involves 
making assumptions about stellar radii and inclinations, 
however, neither of which are well constrained.^^ 

Because of these challenges, our view of the age- 
rotation-activity relation depended until recently on ob- 
servations of handfuls of stars in the field and in a small 
number of well-studied clusters, with the Hyad es be- 
ing a particularly key cluster (e . g., Radick ct al. 119871: 
Jones et al.l Il996t iStauffer et ahl |1997[ iTerndrup et all 



20001 ). Largelv because of the advent of time-domain sur- 



veys, with their emphasis on wide-field, automated, high- 
cadence observing, it is now possible to monitor stellar 
rotation on an entirely new scale (e.g..llrwin et aLll2007l : 
iMeibom et "all 120091 ; [Hartman et alJl2010D. The Palo- 
mar T ransient Factorv fPTF HLaw et al.ll2009l : lRau et all 
120091) provides deep, multi-epoch photometry over a wide 
field-of-view, and our Columbia/Cornell/Caltech PTF 
(CCCP) survey, one of PTF's Key Projects, is leverag- 
ing this capability to measure rotation periods in open 
clusters of different ages. 

Our first CCCP target, Pracsepc,!^ 08 40 24 -1-19 41, 
is a n earby (^180 pc; [ van Lccuwen 2009), rich (~1200 
stars: JKraus fc Hille nbrand 2007), and intermediate-age 
('--^GOO Myr; iDelorme et al.i 20111 cluster that shares 
many characteristics with the Hyades. Until re- 
cently, only five rotation periods — for mid- to late- 
M dwarfs — had been m easured for Praesepe members 
(jScholz fc Eisloffe]|l2007f ) . These periods were often com- 
bined with (less sparse) data for high-mass Hyads in or- 
der to inf er the mass-rotation re lation for 600-Myr-old 
stars (e.g. Jlrwin fc Bouvieill2009f ). This was particularly 
unsatisfying as the Hyades and Praesepe were the two 
oldest clusters with measured rotation periods, and were 
therefore essential in studying the evolution of the age- 
rotation relation from ages of a few 100 Myr to the age 
of the Sun. 

Fortunately , the situation has imp roved significantly in 
the past year. IDelorme et al.l (|201ll ) surveyed the Hyades 
and Praesepe as part of the SuperWASP exoplanet- 
search program. SuperWASP 's sensitivity, tuned to dis- 
cover exoplanets transiting nearby bright stars, enabled 
the measurement of rotation periods for 52 l ate-F to late- 
K/ear ly-M stars in Praesepe. Meanwhile, IScholz et al.l 
(|2011l ) added 49 rotatio n periods ( of which 24 ar e consid- 
ered very robust) to the lScholz fc Eisloffel ()2007l ) sample, 
with the bulk of this new sample being of spectral type 
M3-M5.^^ (Thanks to Kepler, rotation periods have now 
also been mea s ured in a 1-Gyr-old cluster, NGC 6811; 
IMeibom et all (l20ll .) 

We report stellar rotation periods for 40 late-K/early- 

^^ A further limitation of the Doppler broadening technique is 
that it is sensitive only to stars rotating faster than some threshold 
set by the spectral resolution. 

^'^ Also known as the Beehive Cluster and M44. 

^'^ These spectral types are from the'Kraus & Hillenbrand' f2007) 
catalog; a few Scholz & Eisloffel (2007) and Dclormc ct al. (2011) 
stars lack spectral types in this catalog because they are too faint 
or too bright. 



M Praesepe members derived from our first season of 
PTF observations. Our campaign produced ~200 dis- 
tinct observations of four overlapping fields designed to 
inc lude a large number of Prae sepe members identified 
by IKraus fc Hillenbrandl (|2007l ). In Section [3 we de- 
scribe our data and in Section |3| our period-finding al- 
gorithm We also compare our periods with those de - 
rive d bvlScholz fc Eisloffeil (l2007h , IDelorme et all ((20T1 , 
and IScholz et al.l ( 20111 ) for the stars for which they 
also measured Prot, and flag potential binary systems 
among our rotators. In Se ction [H w e comb i ne ou r 
Praesepe data with that of IScholz fc EisloffeJ (|2007[) . 
IDelorme et al.l (|2011| ). and I Scholz et al.. (,2011.1 and com- 
pare color-period relations and mass-period distribu- 
tions derived from these data to those derived from the 

the 150- 



Hyades (using data from IDelorme et al. 2011) , ~j^^^ ^,j^ 
Myr-o ld cl usters M35 and NGC 2516 (|Meibom et all 
(J2009I ) and llrwin et aTJ (2007) , respectively) , and kine- 
matically selected young and old field star populations 
(jKiraga fc StepienI 120071 ) , as wel l as to gyrochrones de- 
rived from the models of lBarnei ()2010D . We conclude in 
Section O The Appendix lists interesting variable field 
stars identified in our Praesepe observations. 

In a forthcoming companion paper we use the results of 
our spectroscopic campaign with the 2.4-m Hiltner tele- 
scope at MDM Observatory and the WIYN 3.5-m tele- 
scope at NOAO, both on Kitt Peak, AZ, to examine the 
relationship between rotation and activity in Praesepe. 
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Fig. 1. — CMD of Praesepe members with 8 < rgDSS < 21 mag. 
Photometric errors are plotted but are typically smaller than the 
symbols. Overplotted in red are members with PTF detections; the 
solid line is the cluster main sequence. Stars with r ^ 14 — 16 and 
(r — K) > 4 are ones for which the SDSS photometry is saturated, 
and for which accurate photometry could not be calculated from 
the wings of the point spread function. 



2. OBSERVATIONS AND DATA REDUCTION 

The Palomar Transient Factory is a transient detection 
system comprised of a wide-field survey camera mounted 
on the automated Samuel Oschin 48 inch telescope at 
Palomar Observatory, CA, (known as the P48), an au- 
tomated real-time data reduction pipeline, a dedicated 
photometric follow-up telescope (the automated Palomar 
60 inch), and an archive of all detected sources. The tech- 
nical aspects of PTF and the project's science goals are 
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TABLE 1 
PTF Observations of Praesepe 



Field 
Number 



Field 
Center 



Number of 
Observations 



110001 08 39 +19 15 185 

110002 08 39 +20 15 189 

110003 08 44 +19 15 195 

110004 08 44 +20 15 198 



descri bed in detail in iLaw et aT] ()2009D and iRau et all 
(|2009f) . The P48 survey camera is based on the CFH12K 
mosaic ca mera formerly at t he Canada-France-Hawaii 
Telescope (jRahmer et al.l2()08l ). The camera has 12 chips 
(one of which is not working), 101 megapixels, 1" sam- 
pling, and a 7.26 deg^ field-of-view. Observations are 
performed in either Mould R or Sloan Digital Sky Sur- 
vey (SDSS) 5, or with a set of Ha filters. Under median 
seeing conditions (1-1") the camera achieves 2.0" full- 
width half-maximum images, and reaches 5 standard de- 
viation (a) limiting AB magnitude s of nin ~ 21.3 and 
mpi Rs 21.0 mag in 60 s exposures (" Law et 311120101 ). As 
of 2011 May, the PTF footprint included 7000 deg^ that 
have been visited at least 25 times, with nearly 1000 deg^ 
having been imaged at least 100 times. 

Four overlapping 3.5 x 2.31 deg fields covering the cen- 
ter of Praesepe were imaged by PTF beginning on 2010 
February 2 and ending on 2010 May 19.^^ Because we 
shared some of our observing time w ith PTF's transiting- 
planet search (see Law et aL|l2"011[ ). there were multiple 
nights early in our campaign and one in April when the 
cluster was observed every 15 min, resulting in 15-30 im- 
ages per night. For most of the campaign, the fields were 
observed 1-2 times a night when the weather allowed, re- 
sulting in close to 200 observations for each field (see 
Table [1} . This observing cadence was sensitive to Prot 
from a few to a few hundred hours, covering the range 
occupied by the fe w cluster members with measured pe- 
riods then known ("Scholz fc Eisloffe]|l2007f ). 

Kraus & Hillenbrand (200^3) combined data from SDSS 
(|York et all 12000). the Two Micron All Sky Sur- 
yev (2MAS S; ISkrutskie eTall [2006 ). and USNO-Bl.O 
(jMonet et al. 2003) to calculate proper motions and pho- 
tometry for several million sources within 7 deg of Prae- 
sepe's center. This census covers a larger area of sky 
and is deeper than any previous proper motion study of 
the cluster. The resulting catalog includes 1129 candi- 
date members with membership probability >50% (here- 
after referred to as the P50 stars); 442 were identi- 
fied as high-probabil i ty ca ndidates for the first time. 
iKraus fc Hillenbrandl (|2007l ) estimated that their survey 
is >90% complete across a wide range of spectral types, 
from FO to Mh}^ 

Of t he 1129 P50 members in the IKraus fc Hille nbrandl 
(|2007f) catalog, 923, or over 80%, he within the CCCP 
footprint. Of these, 661 are fainter than the PTF satu- 
ration limit (--14 mag): PTF detected 534 (or 81%) of 
these candidate members, with the rest falling within 
chip gaps or on the dead chip. Figure [1] is a color- 

^^ The full extent of our observation al footprint was about 18 
deg^; see Figure 1 in lCovev et al.l II2010I ). 

-^40 k nown members are not included in the 
; Hillenb randl 120071) catalog as they are too bright 
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Fig. 2. — a versus median R magnitude of objects detected in 
multiple epochs in field 110003. Praesepe members are shown as 
red stars; the ~30,500 non-members with PTF light-curves in this 
field are shown in black. At the bright end, the scatter in the 
light-curves exceeds the formal photometric errors by factors of a 
few, indicating that the precision is limited by systematic effects 
rather than by random photometric error. We plot in blue the 
median a for non-members obtained when placing the magnitudes 
in bins of width 0.5 mag (e.g., for the 1266 non-members with 
16.5 < -R < 17 mag, the median a = 0.02 mag). Praesepe members 
have systematically higher a (and are thus more variable) than the 
median field star for ij < 17. 
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Fig. 3. — Sample light-curve for the Praesepe member AD 1508 
illustrating the mix of observing cadences. The x-axis is the num- 
ber of Julian days since 2009 January 1; periods of high-cadence 
observing are visible around days 410 and 470. 142 separate pho- 
tometric measurements are plotted; the errors are ~0.02 mag on 
average. 

magnitude diagram (CMD) of Praesepe members de- 
rived from SDSS and 2MASS photometry for stars with 
8 < r < 21 mag. For plotting purposes we apply 
minimal quality cuts, requiring only that that the er- 
rors in__both__r_aiid__K_be^^0.1 mag. 1105 stars in 
the IKraus fc Hillenbrandl (|2007D Praesepe catalog and 
529 of the mem bers detected by PTF meet these crite- 
ria. IKraus fcHill enbrand (2007) provide spectral types 
for Praesepe stars based on spectral-energy-distribution 
(SED) fitting. We apply a polynomial fit to these SEDs 
and thereby define an interpolated single-star main se- 
quence for the cluster. The PTF-detected members are 
late-K through early-M stars, as expected given the dis- 
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tance to Praesepe and the PTF exposure time. 

Aperture photometry was measured for each can- 
didate member at ea ch epoch using SExtractor 
(jBertin fc Arnouta Il996f ). Precision differential pho- 
tometry was per formed using the pipehne described in 
ILaw et all (|2011D . and the base zero-points were defined 
using the SDSS magnitudes of several hundred reference 
stars. The typical long-term photometric stability is ap- 
proximately 5 millimag for the brightest stars, and is 
photon-limited for all stars fainter than i? ~ 15 mag. 
Figure [21 shows the a versus median R for objects de- 
tected in multiple epochs for one of the Praesepe fields. 

Positional matching was used to merge detections 
across epochs, producing a single light-curve for each 
source. Our Praesepe fields overlapped substantially in 
the cluster center; for stars observed in multiple fields, 
light-curves were merged after applying small offsets to 
the light-curves from each field to ensure a comnron me- 
dian magnitude. A sample light-curve for a Praesepe 
member is shown in Figure [31 



correspond to a periodogram power threshold >20. We 
therefore adopted a conservative power threshold of 30 
to select potentially periodic cluster members. 

We then visually inspected the output of our search for 
each candidate. Periodograms were checked to confirm 
the presence of a single narrow peak, well separated from 
the underlying background power; further scrutiny estab- 
lished that the periodic behavior was visible and stable 
across the full light-curve, well-sampled in phase, and of 
an amplitude at least comparable to the observational 
noise. The periodograms and phased light-curves for 
the high-confidence rotators are shown in Figure [TT] and 
listed in Table [H Our analysis produced high-confidence 
measurements of Prot ranging from 0.52 to 35.85 d for a 
total of 40 stars. 37 of these stars have Pm em > 95%, 
as calculated by iKraus fc Hillenbrandl ()2007n , with two 
of the other stars having Pmem > 94%. The remaining 
rotator, JS 634, is a relatively low-probability member of 
the cluster, with Pmem = 62.3%. Radial velocity obser- 
vations will be required to confirm its membership. 



3. PERIOD MEASUREMENTS 

We used a modified version of the Lomb-Scargle al- 
gorithm to search our light-curves for period ic signals. 
We followed Eq. 11 of iFrescura et "all (|2008f) to define 
frequency grids based on the number of measurements 
obtained for each star and the time-span of each light- 
curve. We oversampled by a factor of five to ensure maxi- 
mum sensitivity to periodic variability. Lomb-Scargle pe- 
riodograms were then computed iteratively: an initial pe- 
riodogram was calculated from those data points within 
6(7 of the light-curve's median magnitude. The light- 
curve was phase-folded using the period corresponding 
to the frequency with the maximum periodogram power 
and smoothed with a boxcar window whose width was 
equivalent to 10% of the data points within the full light- 
curve. 

Residuals were calculated for each point with respect to 
this phased, smoothed light-curve. Points with residuals 
greater than 4cr from the smoothed, phased light-curve 
were rejected before calculating a new Lomb-Scargle pe- 
riodogram. After two iterations, the frequency with the 
maximum power was selected as the most likely period 
for that star. Potential beat frequencies between the pri- 
mary periodogram peak and a possible one-day alias, 
typical for ground-based, nigh tly observing campa igns, 
were flagged following Eq. 1 of iMessina et al.l ()2010.1 . In 
this manner we measured potential rotation periods for 
all of the cluster members detected by PTF. 

To test the significance of these period s, we perform ed 
a permutation test on our light-curves (Efron 1982).^^ 
We conducted the analysis described above on each light- 
curve after randomly scrambling the magnitudes mea- 
sured at each epoch. Repeating this test 100 times on 
each scrambled light-curve, we identify the maximum 
measured periodogram peak as the power threshold cor- 
responding to a <1% false alarm probability (FAP) in the 
absence of ordered variations. Across our entire sample, 
this analysis established that a periodogram peak with 
power ^25 corresponded to a FAP <1%; indeed, for only 
three of the 534 stars analyzed here did the 1% FAP 

^'^ For a s imilar approac h, but using randomly generated mag- 
nitudes, see .Frescura et aLl feOOS) . 
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Fig. 4. — CMD with t he Praesepe rotators ident ified by this 
study along with those of IScholz fc Ei sloff'el (2003), IScholz et all 
H2011h . and IDelorme et al.l l|20111) . For cluster stars, we apply 
more stringent photometric tests, e.g., that SATURATED = and 
that the JC-band flag be "A" or "B"; see (Stoughton et al. (2002); 
ISkrutskie et al.l 120061 ). To the single-star main sequence plotted 
in Figure [T] we add a binary main sequence offset by 0.75 mag. 
The inset is a close-up of the region occupied by our rotators; stars 
above the dot-dashed line are candidate binaries (as is JS 497, 
which has an (r — K) = 6). Other cluster stars are omitted here 
for clarity. 



3.1. Comparison to other surveys of Praesepe 

IDelorme et al.l ()2011[ ) surveyed two Praesepe fields 
with the eight cameras of SuperWASP as part of their 
search for transiting exoplanets. This survey produced 
60-70 usable images per night for 60 nights spread over 
130 nights, and these authors searched for rotation pe- 
riods between 1.1 to 20 d for cluster members within 10 
deg of the cluster center. As SuperWASP is optimized to 
detect transits around nearby, bright stars, the 52 stars 
for which this survey measured rotation periods (46 of 
which have Pmem > 95%) are significantl y brighter than 
those in our sample, and there is only one IDelorme et aLl 
(|2011. ) star, JS 545, for which we also measured Prot- 
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TABLE 2 
High confidence rotators 



Name 



RA 



Dec 



SpT 



Mf 



Ml 
(Mq) 



M2 



Ave. PTF 
R (mag) 



#of 
Obs. 



(d) Power 



HSHJ 428 


08 42 37.61 


+19 59 18.91 


M3.5 


99.5 


7.10 


0.27 


0.31 


17.90+0.05 


358 


1.69 


69.60 


JS 687'' 


08 35 59.44 


+20 04 40.54 


M3.7 


99.8 


6.71 


0.33 


0.37 


17.55+0.03 


172 


1.76 


42.33 


JS 718=^ 


08 40 04.16 


+19 24 50.27 


M3.4 


99.7 


6.68 


0.33 


0.37 


17.37+0.03 


448 


2.62 


113.45 


JS 667=^ 


08 33 17.98 


+19 16 32.75 


M3.7 


98.9 


6.52 


0.36 


0.40 


17.40+0.03 


347 


2.30 


102.98 


HSHJ 87=^ 


08 35 47.22 


+18 08 29.97 


M3.2 


98.8 


6.41 


0.38 


0.41 


17.04+0.02 


168 


2.12 


41.65 


2MASS 08505688+1936579 


08 50 56.87 


+19 36 57.89 


M2.5 


97.1 


6.40 


0.38 


0.41 


16.76+0.02 


320 


11.95 


49.70 


AD 3814'^ 


08 50 49.84 


+19 48 36.51 


M3.4 


97.9 


6.36 


0.39 


0.42 


17.01+0.02 


363 


7.43 


65.47 


AD 2552 


08 39 22.44 


+20 04 54.77 


Ml. 9 


99.9 


6.22 


0.41 


0.44 


16.33+0.01 


530 


25.36 


127.43 


JS 241'^ 


08 38 37.47 


+19 15 28.67 


M3.2 


99.8 


6.11 


0.43 


0.46 


16.66+0.02 


526 


0.52 


55.32 


JS 365 


08 40 30.58 


+19 55 58.80 


M2.7 


99.8 


6.10 


0.43 


0.46 


16.42+0.01 


537 


17.04 


153.22 


JS 644 


08 47 09.10 


+18 11 37.30 


M2.7 


97.3 


6.13 


0.43 


0.46 


16.37+0.01 


180 


8.32 


40.63 


JS 95=^ 


08 35 40.15 


+18 42 28.32 


M3.3 


99.4 


6.07 


0.44 


0.46 


16.74+0.02 


168 


1.98 


46.04 


JS 250=^ 


08 38 51.03 


+19 51 02.03 


M3.0 


99.5 


6.07 


0.44 


0.46 


16.60+0.02 


534 


35.85 


75.85 


JS 525 


08 43 12.92 


+18 31 50.83 


Ml. 9 


99.2 


5.99 


0.45 


0.48 


16.18+0.01 


178 


22.60 


38.89 


JS 506=^ 


08 42 52.28 


+19 51 45.92 


M2.5 


99.4 


5.92 


0.47 


0.49 


16.27+0.01 


364 


3.97 


101.18 


JS 283 


08 39 16.79 


+19 47 42.63 


Ml. 4 


99.8 


5.79 


0.49 


0.51 


15.71+0.01 


534 


19.62 


116.87 


JS 110 


08 36 08.97 


+19 13 48.02 


Ml. 3 


99.7 


5.75 


0.50 


0.51 


15.68+0.01 


170 


18.21 


40.71 


JS 237 


08 38 32.83 


+19 46 25.61 


Ml. 2 


99.9 


5.72 


0.51 


0.52 


15.59+0.01 


535 


18.57 


88.75 


JS 148=^ 


08 36 48.95 


+19 18 59.30 


Ml. 8 


99.8 


5.66 


0.52 


0.53 


15.78+0.01 


532 


1.17 


189.58 


JS 415 


08 41 10.52 


+18 16 07.05 


Ml. 7 


99.3 


5.67 


0.52 


0.52 


15.75+0.01 


356 


9.20 


89.30 


JS 457 


08 41 50.05 


+19 39 34.68 


Ml. 4 


99.7 


5.60 


0.53 


0.53 


15.43+0.01 


532 


18.90 


140.76 


JS 468 


08 41 59.35 


+19 44 45.11 


Ml. 2 


99.8 


5.60 


0.53 


0.53 


15.42+0.01 


530 


16.77 


73.70 


JS 174=^ 


08 37 19.91 


+19 03 11.92 


M3.0 


99.7 


5.54 


0.54 


0.54 


15.80+0.01 


347 


2.03 


61.36 


JS 187 


08 37 32.42 


+19 31 17.98 


Ml. 9 


99.0 


5.50 


0.55 


0.55 


15.42+0.01 


372 


1.81 


113.91 


JS 46 


08 33 50.76 


+19 46 58.62 


M0.5 


99.5 


5.44 


0.57 


0.56 


14.98+0.00 


347 


17.00 


80.65 


JS 364 


08 40 28.23 


+18 56 08.96 


K7.8 


99.2 


5.42 


0.57 


0.56 


15.08+0.00 


356 


17.33 


60.84 


JS 489 


08 42 34.86 


+20 59 40.80 


Ml.O 


95.9 


5.46 


0.56 


0.56 


15.18+0.00 


183 


17.33 


61.77 


JS 505 


08 42 49.68 


+18 51 35.10 


Ml. 5 


94.4 


5.45 


0.56 


0.56 


15.33+0.00 


176 


13.50 


49.98 


JS 430=^ 


08 41 24.46 


+20 07 49.52 


M2.5 


99.6 


5.35 


0.58 


0.57 


15.78+0.01 


538 


1.40 


95.15 


JS 634'''b 


08 46 38.21 


+19 52 44.78 


M2.5 


62.3 


5.22 


0.61 


0.59 


15.42+0.01 


359 


3.24 


80.98 


JS 181 


08 37 28.45 


+20 36 28.54 


M0.8 


99.1 


5.20 


0.62 


0.60 


14.84+0.00 


185 


7.22 


82.29 


JS 432 


08 41 24.18 


+18 14 02.69 


Ml.O 


99.3 


5.17 


0.62 


0.60 


15.02+0.00 


356 


2.09 


134.68 


JS 545 


08 43 44.73 


+21 12 34.32 


M0.8 


97.5 


5.12 


0.63 


0.61 


14.60+0.00 


183 


2.83 


53.99 


2MASS 08492676+1831196 


08 49 26.76 


+18 31 19.54 


MO. 5 


97.8 


5.10 


0.64 


0.61 


14.65+0.00 


178 


9.20 


64.69 


KW 563'''<= 


08 39 09.86 


+19 46 58.93 


M2.0 


99.8 


5.06 


0.65 


0.62 


15.30+0.01 


539 


4.33 


132.19 


KW 570 


08 39 41.03 


+19 59 28.81 


M0.4 


99.9 


5.03 


0.65 


0.62 


14.59+0.00 


465 


4.18 


189.34 


JS 497=^ 


08 42 42.07 


+19 17 32.26 


K7.3 


94.4 


4.76 


0.71 


0.67 


14.09+0.01 


153 


12.68 


50.40 


HSHJ 15=^ 


08 31 40.45 


+19 47 54.19 


M2.0 


95.2 


4.74 


0.71 


0.68 


14.36+0.00 


347 


9.36 


114.60 


JS 159=^ 


08 36 56.25 


+18 57 47.97 


M0.9 


99.6 


4.62 


0.74 


0.70 


14.36+0.00 


275 


4.81 


95.55 


AD 1508=^ 


08 31 29.87 


+20 24 37.49 


MO.l 


98.3 


4.48 


0.77 


0.74 


14.07+0.00 


142 


1.55 


63.48 



























vey resulted in rotation periods for 49 stars, all but one 



of which have cluster P„ 



> 95%. These stars are 



Note. — Spectral types and membership probabilities are taken from iKraus &: Hillenbrandl J2007D : for the former, the formal uneertainty is 
0.1 spectral classes. However, the systematic uncertainty in the underlying definition of spectral types is ~0.5 spectral classes for M dwarfs, and 
this systematic uncertainty will be reflected in the color-spec tral type relations us ed for SED fits. Mk is calculated assuming a distance of 181.5 
pc (van Leeuwen 2009). Mi is esti mated using th e empi rical IDelfosse et all ll2000f) relation, while M2 is estimated using the theoretical model of 
iDo tter et al. (2008). Altho ugh the IDelfoss e et al.l (1200011 relation extends to stars with M^^ — 4.5 mag, the predicted mass values diverge by up 
to about 5% from those of IDotter et aU ( 120081) for stars brighter than Mk — 5.5 mag. The number in bold is the mass we assigned for a given 
star. We provide the mean magnitude of each light-curve after filtering on flags has been applied; the quoted error is the mean of the Icr errors 
on the magnitudes used in calculating this mean magnitude. 
^ Candidate binary system. 

b P„em = 62.3%; all other stars have Pmsm > 94% (IKraus fc Hillenbrandll2007l) . 
'^ Listed as the variable V*AXCnc in SIMBAD. 

The overlap betwee n o ur sample and t hose of 
iScholz fc Eisl5ffeil (l2 007l) andlSc holz eFaLl (|20ll is more 
significant. iScholzfc Eisloffcl (2 007f ) published the first 
rotation periods for Praesepe, for five mid- to late-M- 
type cluster members. These data were based on 125 
images of a 0.36 deg^ field taken over the course of about 
a month in early 2001 with the 2-m Schmidt telescope at 
the Thiiringer Landessternwarte Tautenburg (Germany), 
and on 108 images obtained over the course of a week two 
years later with the 3.5-m telescope at Calar Alto Ob- 
servatory (Spain) of a different 0.25 deg^ Praesepe field. 
One star for which they measured Prot, JS 687, appears 
i n our sample o f rota tors. 

IScholz et aU (|20lH ) observed four 0.32 deg^ Praesepe 
fields around the cluster center with the 2.5-m Isaac New- 
ton Telescope (Spain), obtaining a few hundred images 
for each field over the course of eight nights. This sur- 



typically fainter than those in our sample, but there are 
seven for which we measured Prot- 

The nine stars with published Prot values for which we 
have measured a rotation period are listed in Table [H 
The agreement between our values for these stars' pe- 
riods and those in the literature is excellent, with dif- 
ferences >,2% only for the two stars with Prot > 4 d, 
KW 563 and KW 570. The se are among the s tars with 
the longest Prot reported bv IScholz et al.l (|2011 f): neither 
is among the 24 stars flagged as having the most robust 
Prot measured by these authors, likely because of the dif- 
ficulty of accurately measuring periods longer than half 
the duration of a monitoring campaign. 

Our data occupy a unique area in mass-period space. 
The stars for which we measure Prot ai'e later than those 
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TABLE 3 

Stars with previously reported periods 



Name 



SpT 



M (MCd) 
Lit. PTF 



Proi 

Lit. 



(d) 



PTF 



HSHJ 428 M3.5 0.26 0.27 



JS 687^ 
JS 718=^ 
JS 241=^ 
JS 506=^ 
JS 430=^ 
JS 545 
KW 563^ 
KW570 



M3.7 
M3.4 
M3.2 
M2.5 
M2.5 
M0.8 
M2.0 
M0.4 



0.29 
0.33 
0.42 
0.47 
0.54 

0.60 
0.60 



0.33 
0.33 
0.43 
0.47 
0.57 
0.61 
0.62 
0.62 



1.71 
1.76 
2.65 
0.53 
3.97 
1.37 
2.83 {0.004)*= 
4.85 
4.27 



1.69 
1.76 
2.62 
0.52 
3.97 
1.40 
2.83 
4.33 
4.18 



Note. 



Data for JS 687 are from 



Eisloffell 



m I SchoIz c- . 

I I2007I) and fo r JS 545 frorn IDelorme et al'l I I2011I) : all other 
data are from IScholz et al.l (120111) . 
^ Candidate binary system. 

P^ot is the average of the two periods whose difference is 
indicated in parentheses. 

identified as rotators by IDelorme et al.l (120111) and gen- 
erally earli er than those identi fied by IScholz fc EislofFel 
((200'7i) and IScholz et all (|2011h . fihing the gap between 
solar-type and mid-to-late M dwarf rotators in Praesepe. 
Furthermore, the cadence and time-span of our observa- 
tions gives us sensitivity to the fast rotators identified by 
both of these groups an d to the slow rotators identified 
at the high-mass end bv IDelorme et al.l (|2011h . 

3.2. Identifying potential binary systems 

The presence of binary systems could affect our in- 
terpretation of the Praesepe mass-period relation in two 
ways. First, treating a binary system as though it is a 
single star will lead to an erroneous mass estimate. Sec- 
ond, close binaries might be tidally locked and rotating 
faster than single stars of the same spectral type.^* 

We use the cluster CMD to identify candidate binary 
systems among our rotators. iSteele fc JamesonI ()1995D 
showed that in the Pleiades the effect of binaries on the 
CMD is to create a second "main sequence" lying above 
that of single stars and offset by ~0.75 mag for a given 
color; we apply this same o ffset to the main sequence we 
derive from our fit to the iKraus fc Hillenbrandl ()2007f ) 
stellar SEDs (see Figure S]). Stars along this sequence are 
likely to be in close to equal mass systems, with systems 
with m ore extreme mass ratios lying in the interven ing 
region (ISteele fc Jamesonlll995l: fHodgkin et al.lll9990 . 

We adopt the method of iHodgkin et al.l ()1999[ ) and 
identify 18 candidate binary systems lying above the mid- 
point between the single-star and binary main sequences. 
These are highlighted in Figure [H and in Table [3 We 
note that the corresponding binary fraction is consistent 
(if slightly belo w) what has been r eported previously for 
the cluster (e.g. . IScholz et al.ll2011h . Radial velocity mea- 
surements are required to confirm that these are binaries. 

4. DISCUSSION 

4.1. Color-period relation for solar-type stars & 
comparison to younger clusters 

In Figure [51 we show the color-period distribution for 
Praesepe rotators, adding our new PTF data to the data 

^® Some spot configurations coul d also lead to errors in our Prot 
measurements. As pointed out bv IScholz et al.l (120111 1. two spots 
of the same size separated by 180 deg in longitude would lead to a 
measurement of half the true stellar rotation period, for example. 
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Fig. 5. — Color-p eriod distribution for Praesepe, with the 
IDelorme et al.l 120111) color-period relation overplotted. 
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Fig. 6.— Color-per iod distribution for M35 and NGC 2516. The 
IDelorme et al.l 120111 ) color-period relation shown in Figure \E\ has 
been spun up by (600/150)'''^ , following Skumanic^ 119721) . to pro- 
duce the black color- period relation. The red line takes the same 
IDelorme et"all 120111) relation and spins it up by (600/150)'''^^. 



collect ed bv IScholz fc Eisloffell (120071 ). IDelorme et all 
(|2011f) . and IScholz et al.l (|2011D . These data confirm 
that the color-period relation for Praesepe can be consid- 
ered sing le- valued for stars blu er than {J — K) ^ 0.8, as 
found bv IDelorme et all (|2011|) . For redder, lower- mass 
stars, however, the distribution of stars in color-period 
space is scattered, with populations of fast and slow ro- 
tators. In particular, the PTF data allow us to measure 
Prot > 10 d in stars with 0.8 < (J - K) < 0.95, so that 
the sharp break in the color-period distribution noted by 
IScholz et al.l (120111) is no longer evident. ^^ 

[Scholz et al.' (2011) used the data collected by 
EWin et al. (2007) for the 150-Myr-cluster NGC 2516 to 
test models of rotational evolution by comparing the dis- 
tribution of rotational periods in their Praesepe data to 
that for the younger cluster. We sought to replicate this 
analysis for bot h high- and low-m ass stars; to this end 
we added to the llrwin et al.l ()2007f ) data the rotation pe- 

^^ Long-term monitoring of the faintest Praesepe members is 
needed to determine whether the absence of slow rotators with 
(J — ii')>1.0is real or an observational bias. 
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riods collected bv lMeibom et al.l (|2009f) for M35, another 
'--^ISO-Myr-cluster. For both clusters we required that the 
2MASS photometry have either an "A" or "B" quality 
flag; this resulted in a sample of 300 stars for M35 and of 
294 stars for NGC 2516. The resulting color-period plot 
for the two clust ers is showii i n Fig ure El 

We tested the iSkumaniciil (|1972l ) law for the higher- 
mass stars in both clust ers by spinning up the color- 
period relation derived bv lDelorme et al] (|2011f ) for Prae- 
sepe by (600/150)"'^. The resulting color-period relation 
is plotted in Figure [5] This spun-up color-period rela- 
tion is fairly insensitive to the exact ages used for the 
clusters, and does not appear to describe the M35 stars 
well regardless of what (reasonable) age combination is 
used for the two clusters. Our analysis suggests that the 
the age dependence is closer to t^-^^ for the slow rotators 
with 0.5 < {J — K) < 0.8; furthermore, it is difficult to 
find a single-valued color-period relation that describes 
these slow rotators and their M35 neighbors bluer than 
(J-ivr)~0.5. 

While the spin-down derived from this co mparison is 
less st eep tha n that found by apply i ng the iSkumanichI 
(|1972l) law, iCollier Cameron et all (|2009D find that 
Prot oc i°'^^ for (slightly earlier) solar-mass stars between 
the ages of the Hyades/Coma Ber (another ~600-Myr- 
old cluster) and the age of the Sun. This implies that 
these stars are spinning down faster once they reach 
ages >600 Myr, and is con sistent with e.g., the t heo- 
retical picture presented bv llrwin fc Bouvien (|2009( ): if 
the radiative core rotates more rapidly than the convec- 
tive envelope when the star is very young, the transfer of 
angular momentum from core to envelope will "soften" 
the rotational braking of the surface layers early on. This 
will produce a sub- Skumanichj (jl972l ) spin-down between 
^--^150 and 600 Myr, with the transition to a steeper spin- 
down rate occurring only once the transfer of angula r 
momentum ends (see Figure 5 of llrwin fc Bouvierll2009f ) . 

For the color range we were most sensitive to with 
our PTF observations, however, working in color-period 
space was somewhat unsatisfying, as our sample spans a 
rel atively narrow range of ( J — K) colors (as d o tho se 
of IScholz fc Eisl5ffeil pOOl and IScholz et al.1 (fMl V 
We therefore chose to examine the relationship between 
masses and periods for these samples; this was further 
motivated by the fact that mass, unlike color, does not 
change significantly as the star evolves from the pre-main 
sequence onward. 

4.2. Masses, gyrochrones, & mass-period comparisons 

A number of previous studies have used a star's (B—V) 
as a proxy for its mass, and by calibrating the resulting 
color-period sequence in clusters of various ages, derived 
analytic expressio ns for a star's Pmt as a function of (B — 
y) a nd age (e.g., iBarnesI 120031 : iMamaiek fc Hillenbrandl 
|2008[ ). Many of the lowest mass Praesepe members lack 
accurate (B — V) colors, however. Furthermore, placing 
the Praesepe rotators and rotators from other clusters 
on a common mass scale requires a mass proxy with a 
large dynamic range. As shown in Figure [H Prot has 
been measured for Praesepe members spanning ~13 mag 
in r, an achievement only made possible by combining 
light-curves obtained by facilities with distinctly different 
sensitivities. 



Fortunately, this span is somewhat compressed in the 
near infrared, as demonstrated by the range of (r — K) 
colors covered by the Praesepe cluster sequence. To ex- 
tend our analysis to the low-mass Praesepe cluster mem- 
bers, we therefore use mass estimates based on 2MASS 
photometry, which is available for all of the clusters 
we co nsider. Furthermore, as noted by iDelfosse et al.l 
(|2000l ). near- infrared absolute magnitudes are better cor- 
related with mass than their optical equivalents, at least 
for M dwarfs. 

We began by calculating each star's absolute K mag- 
nitude (My), using a distance to Praesepe of 181.5 pc 
(|van Leeuwen|[2009r i. The empirical M^^-mass relation is 
the best calibrated of the absolute magn itude- mass re- 
lations described bv IDelfosse et al.l ()2000[ ). and for stars 
with Mk > 5.5, we derived masses using this relation. 
For stars with Mk < 5.5, we estimated masse s using the 
theoretical M/^-mass relation developed by Potter et al.l 
(2008) for a 600 Myr, solar-metallicity population. This 
relatio n agrees well with the em pirical M/^-mass rela- 
tion of iHenrv &: McCarthy! (|1993D but extends to higher 
masses. For completeness, both mass estimates are listed 
in Table [2] for each of our rotators; the adopted mass es- 
timate is highlighted in boldface. 

We followed the same steps to estimate the masses 
of Praesepe members with r ot ation periods m e asure d 
by IS cholz fc EislHlea (l2007h . IDelorme et all ((20T1 . 
and Scho lz et al.l ()2011[ ). Comparing our mass esti- 
mates to tho se der ived by Scholz & Eisloffcl (2007) and 
IScholz et aP ()2011[ ) from isochrone models indicates that 
the two techniques produce mass estimates that are 
consistent to within a few percent. Errors in the as- 
sumed cluster distance propagate linearly into the de- 
rived masses, so the uncertainties in our mass esti- 
mates reflect the ~5% uncertainty in the cluster distance 
jvan Lccuwcn 2009); systematic uncertainties in the 
[Pelfosse et al.. (,2000. ) relation are also of order ~5— 10%, 
and we therefore adopt 10% as the typical uncertainty 
in our derived masses. The location of the full sample in 
mass-period space is shown in Figure [3 

4.2.1. Comparisons to model gyrochrones 

A major motivation for surveying stellar rotation in 
open clusters is to calibrate the relationship between a 
star's rotation period and its age. Wc have therefore 
used the formalism of 'Barne s fc KinJ (J2010) and Earned 
(2010) to compute gyrochrones that quantify this rela- 
tionship for a given age and zero age main sequ ence 
(ZAMS) rotation period. Eq. 22 in IBarnesI (|2010f ) re- 
lates a star's convective turnover time (r) to a quadratic 
function of the star's age and the ratio between its ZAMS 
and present day rotation pe riod. Adopting an age of 600 
Myr for Praesepe, we follow iBarnesI (|2010f ) in calculating 
gyrochrones for stars with ZAMS periods ranging from 
3.4 days down to 0.12 days, corresponding to break-up 
for a solar- type star. 

Wh en age and ZAMS period are fixed, the IBarnesI 
(2010) quadratic equation becomes a double-valued re- 
lationship between r and current-day period. Using 
the relation between stellar ma ss and global convect ive 
turnover timescale tabulated bv lBarnes fc Kirn! ()2010f ) to 
project the r-period relationships into mass-period space, 
we overlay in Figure|7]the resulting 600 Myr gyrochrones. 
These gyrochrones bound the region of the mass-period 
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Fig. 7. — Mass-period distribution for Praescpe, with model gy- 
rochrones from Barnes (2010) ovorplotted. Each gyrochrone corre- 
sponds to the predicted position at 600 Myr of stars of a range of 
masses but with the same ZAMS rotation period, Po- 
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Fig. 8. — Mass-period distribution for the Hyades, again with 
gyrochrones for a 600 Myr population overplotted. 

plane containing Praesepe rotators reasonably well, indi- 
cating a tight mass-period relation for masses >0.7 Mq 
and a significantly broader range of periods at masses 
<0.6 Mq. The sample of Praesepe rotators, however, 
includes several stars that lie outside the region of the 
mass-period plane enclosed by these gyrochrones, sug- 
gesting either that Praesepe stars had a broader range 
of ZAMS Prot than 0.12-3.4 days, or that this analytic 
description does not yet fully replicate the morphology 
of the ^600 Myr mass-period relation. 

4.2.2. Comparison to the Hyades 

We have collected observational data for candidate 
members of the Hya des whose Prot were measured by 
iDelorme et alj (|2011[ ). We estimated masses for these 
stars following the technique described a bove and us- 
ing a distance to the Hyades of 46.45 pc (jvan LeeuwenI 
|2009( ). The Hyades is sufficiently close that many mem- 
bers have direct parallax measurements: for those, we 
adopted the distanc es implied by their in dividual paral- 
laxes as reported bv lDelorme et al.l ()2011[ ). The resulting 
mass-period plot is shown in Figure [51 

Figures[7]and[5]suggest that the transition from a slow- 



rotating to mixed fast- and slow-rotating population oc- 
curs at a lower mass in the Hyades than in Praesepe. 
While the data for low-mass stars are sparser than for 
Praesepe, there are comparatively very few fast rotators 
in the Hyades above a mass of ^0.5 Mq. This strength- 
ens the iDelorme et al.l (|2011) conclusion that the clus- 
ters differ in age by ^^50 Myr, with Praesepe being the 
younger of the two. 

4.3. The evolution of the mass-period relation 

These comparisons of a stellar population to gy- 
rochrones generated for a given age are one test of semi- 
empirical models' ability to describe the angular momen- 
tum evolution of stars. Another is to compare the model 
predictions for a star of a given mass to the data avail- 
able for stars of that mass at various ages. We selected 
Praesepe stars not flagged as potential binaries, sepa- 
rated them into two mass bins, 0.3 < M < 0.5 M ^ and 
0.5 < M < 0.7 Mq, and, following IScholzet all ([2011), 
calculated the median and lO"^ and 90'*^ percentile Prot 
for these two samples {Prot = 3.98, 0.84, and 19.62 d, re- 
spectively, for the low-mass bin; Prot = 13.09, 1. 81, and 
18.21 d for the high-mass one). We then used the lBarn"ei 
()2010f ) models to flnd the corresponding ZAMS Po for 
each of these representative 600 Myr stars, to which we 
assign masses of 0.4 and 0.6 Mq. This Po was fed back 
into the models to predict the Prot of these representa- 
tive stars at ages ranging from 125 Myr to 10 Gyr, and 
including 150 Myr, 650 Myr, 1.5 Gyr, and 8.5 Gyr. 



100.0 E 



10.0 



1.0 



0.1 



M35 
JGC 2516 


Hyades young disk 
Praesepe 


_.-♦-■ 


4- ;'< 


• 


*• y' ' 



I 

old disk 



. - P. = 0.15 d, Barnes 2010 
... P^ = 0.73 d 
■■- P. = 8.61 d 
0.3 < M < 0.5 M. 



0.1 



1.0 
Age (Gyr) 



10.0 
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different ages. The red star is the median Praesepe rotator for 0.3 
< M < 0.5 Mq, while the blue stars are the median and 10*'' and 
go**" percentile rotators. The Po for these stars and evolutionary 
tracks are calculated usin g the models of Barnes (2010). Candidate 
binaries flagged in S I3.2l are not included in these calculations but 
are plotted as empty symbols. 

The resulting evolutionary tracks are plotted in Fig- 
ures ini and [ini along with period data for these mass 
bins from M35, NGC 2516, the H yades, and young 
and old disk stars (1.5 and 8.5 Gyr; iKiraga &: St^piehl 
|2007[ ). The mass estimates for the 150-Myr-old clusters 
were done in the manner described in § 14.21 ^° For the 

20 llrwin et all I I2007I) provide mass estimates for NGC 2516 
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Fig. 10. — Same as in Figure[9] for 0.6 Mq stars in Praesepe. 

iKiraga h Stepieril (|2007[ ) stars, we used the mass esti- 
mates these authors provide , which are based on the em- 
pirical [Delfoiii^er^ (|2000D relation for My. 

The agreement between the models for a 0.4 Mq star 
and the observed distribution of rotation periods is good 
between the ages of 150 and 600 Myr. The periods pre- 
dicted for the Praesepe median and 10*^ and 90**^ per- 
centile Prot stars at 150 Myr fall within the range of pe- 
riods observed in M35 and NGC 2516 for stars between 
0.3 and 0.5 Mq, although they are close to the extreme 
period s observed in t hese clusters. 

The iBarnesI ()2010( ) models predict a faster spin-down 
as stars reach an age >1 Gyr than what is seen: by 1.5 
Gyr, the models have converged significantly, while the 
spread in the observed Prot remains roughly equivalent 
to that seen at 600 Myr. Extending the models to the 
age of the thick disk, the rotation periods predicted for 
all percentiles are a few times the (few) periods measured 
for old disk stars. However, the Kira ga fc Stepiehl ()2007l) 
sample is selected from known X-ray-luminous stars, and 
this may bias this sample toward faster rotators. In ad- 
dition, the young and old disk samples are not single-age 
populations; their age spreads are such that they likely 
include stars spanning the full range of ages shown here. 
If the fast rotators are systematically younger than the 
slow ones, the true disagreement with the models is not 
as significant as suggested by Figured 

The agreement between the models and the observed 
spread in Prot for the representative 0.6 Mq stars is also 
good between 150 and 600 Myr, with one clear excep- 
tion. Replicating the 10*^ percentile rotator in Prae- 
sepe requires a ZAMS 0.6 Mq star with Po = 0.05 d, 
which implies that this star is rotating faster than its 
breakup P ~ 0.06 d (assuming that a 0.6 Mq star has R 
= 0.56 Rq\ iDotter et al.ll2008f ). At later ages, the mod- 
els converge even more than for the less massive star, 
so that the match between predicted period distribution 
and that seen at ages >1 Gyr is worse. However, here 
again, the limitations of the disk star sample make com- 
parisons to the models difficult to interpret. 

5. CONCLUSION 



Using PTF observations of four overlapping fields, we 
have measured rotation periods for 40 late-K to mid-M 
stars belonging to the nearby, rich, ~600 Myr open clus- 
ter Praesepe. Our data occupy a unique area in mass- 
period space: we measur e Pmt for stars later t han those 
identified as rotators bv lDelorme et a1 (120111) and gen- 



based on comparisons to the models of lBaraffe et al.l (|1998I 1 . These 
two sets of estimated masses agree to within a few hundredths of 
a solar mass, with our masses being typically slightly lower. 



erally earli er than those identi fied by iScholz fc Eisloffel 
(|200'7[) and IScholz et all (|2011h . filling the gap between 
solar-type and mid-to-late M dwarf rotators in Praesepe. 
Furthermore, the cadence and time-span of our observa- 
tions gives us sensitivity to the fast rotators identified by 
these groups and to the slow ro t ators identified at the 
high-mass end bv lDelorme et al] ()201ir ). 

Our measurements indicate that Praesepe's mass- 
period relation undergoes a transition from a well-defined 
singular relation to a more scattered distribution of both 
fast and slow-rotators at masses ^0.6 Mq, corresponding 
roughly to a spectral type of Ml. The location of this 
transition is broadly consistent with expectations based 
on observations of younger clusters and the assumption 
that stellar-spin down is the dominant mechanism influ- 
encing angular momentum evolution at ^600 Myr. 

A comparis o n to the data recently published by 
iDelorme et al.l (|2011[ ) for the Hyades, widely assumed 
to be coeval to Praesepe, suggests that this transition 
occurs at different characteristic masses in the two clus- 
ters, providing further evidence that Praesepe is the 
younger of the two c luster s. Furthermor e, by using the 
IBarnes fc Kiinl ("2010) and Barnes! ([2010) formalisms to 
evolve the Praesepe Prot in time and comparing the pre- 
dicted Prot with measured Prot in M35 and NGC 2516 
(~150 Myr) and for young and old field star populations 
(1.5 and 8.5 Gyr), we find that stellar spin-down may 
progress more slowly than described by these relations. 

The fixed age mass-period relation is but one projec- 
tion of the underlying stellar age-rotation-activity rela- 
tionship. Previous studies of stellar activity, in clus- 
ters and the field, have derived relationships between a 
star's age and observational tracers of its coronal or chro- 
mospheric activity ( e.g., Skumanich 1972; Radick et al.l 
Il987t iSoderblom et all 12001( 1. Recent studies of chro- 
mospheric activity in low-mass field stars have inferred 
activity lifetimes as a function of spectral type by mod- 
eling the vertical gradient in Ha emission strengths as 
a con sequence of dyna mical heating in the Galactic disk 
(e.g.. lWest et al.ll2008l ). These studies predict that stars 
with spectral types of M2 or later have activity lifetimes 
>1 Gyr. The activity lifetimes of MO-Ml stars are some- 
what less well known, as few active early M stars are 
observed in the field, but appear to be <600 Myr. 

These relations would thus predict that the bound- 
ary between Ha active and inactive Praesepe members 
should occur in the MO/Ml spectral range. The agree- 
ment between the implied mass of the active/inactive 
boundary in Praesepe, near 0.6 M©, and the similar char- 
acteristic mass for the transition from a singular mass- 
period relation to a more scattered distribution of rapid 
and slowly rotators, strengthens the case for an underly- 
ing rotation-activity relation in this cluster. In a forth- 
coming paper, we use the results of our spectroscopic 
campaign with the 2.4-m Hiltner telescope at MDM Ob- 
servatory and the WIYN 3.5-m telescope at NOAO, both 
on Kitt Peak, AZ, to examine this relationship between 
rotation and activity in Praesepe. 
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TABLE 4 
Other Interesting Variable Stars 






Ave. PTF 


#of 


Prot 






SDSS J 


R (mag) 


Obs. 


(d) 


Power 


Type 


083113.96+194951.5 


15.99 + 0.01 


327 


0.26 


136.59 


Eclipsing binary? 


083238.89+210424.6=" 


17.44 + 0.03 


173 


0.55 


57.20 


RR Lyrae 


083426.37+202040.9 


15.53 + 0.01 


348 


0.13 


60.69 


W Uma? 


083525.00+194659.7 


15.96 + 0.01 


170 


0.80 


63.81 


RR Lyrae 


083706.80+185556.2 


16.14 + 0.01 


348 


0.30 


145.78 


W Uma? 


083816.78+182724.6=" 


17.36 ± 0.03 


348 


0.29 


108.32 


RR Lyrae 


084041.70+201612.1 


16.33 + 0.02 


377 


0.13 


123.29 


WUma? 


084214.42+195819.7 


16.33 + 0.01 


531 


0.25 


197.15 


Eclipsing binary? 


084803.57+191340.5 


16.21 + 0.01 


180 


0.49 


55.37 


Known RR Lyrae*" 


084920.49+182617.3 


15.85 + 0.01 


181 


0.53 


63.28 


Known RR Lyrae'"-'" 


084950.50+193232.0 


14.70 + 0.00 


337 


0.18 


137.16 


WUma? 


085112.62+184344.3 


15.06 + 0.00 


181 


0.19 


84.19 


WUma? 



^ Spectrum is available from SDSS. 
^ Classified as RR Lyrae in SIMBAD. 
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performed a search for high-confidence variables within 
this full catalog of light-curves. Candidate variables were 
identified by computing the ratios of the a of the raw 
light-curve to the a of light-curves that were boxcar 
smoothed over windows spanning 9 and 36 epochs. Box- 
car smoothing light-curves dominated by random noise 
will reduce random fluctuations and shrink the a by ^yn. 
By contrast, light-curves dominated by structured vari- 
ability will not converge with y/n, so that the ratio of the 
as of smoothed and raw light-curves serves as a simple 
tool for identifying structured variability. As the light- 
curves presented here are irregularly sampled, the boxcar 
windows corresponded to different timescales at different 
points in the light-curve, and were sensitive to variability 
over a larger range of timescales than if the light-curves 
were regularly sampled. 

By computing a ratios for stars with i? < 18 mag, we 
identified ^6000 stars as candidate variables. We pro- 
duced periodograms in the same manner as described 
in § [3] and defined stars with periodogram peaks >50 
as are likely periodic variables. Visual inspection then 
identified 12 robust detections of large amplitude (>0.1 
mag) variable stars. Phased light-curves for these stars 
are presented in Figure [12] and the stars are tabulated 
in Table S) These are nearly all short period (P < 1 
d) variables; classifying these on the basis of their light- 
curve shapes and amplitudes, we identify five as likely RR 
Lyrae stars (two of these have previously been identified 
as RR Lyrae) and two as candidate eclipsing systems. 
The remaining five show large amplitude sinusoidal vari- 
ations, and on the basis of their short periods, we ten- 
tatively classify these objects as W Uma systems, but 
further study is necessary to confirm their status. 



APPENDIX: INTERESTING VARIABLE STARS IN THE 
PRAESEPE FIELDS 

Our light-curve analysis was largely restricted to stars 
identified p reviously as candidate Praesepe members by 
iKraus fcHil lcnbrand (2007). Candidate Praesepe mem- 
bers total less than 1% of the objects in our target fields, 
however, so that many other variable stars are likely to 
be present in the full catalog of PTF light-curves. 

To explore the population of variable stars that can 
be detected by PTF at moderate Galactic latitude, we 
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Fig. 12.— Other interesting periodic variables in the Praesepc fields. SDSS J084803.57+191340.5 and 084920.49+182617.3 arc both 
previously identified RR Lyrae stars. 
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